Further study of neutrino oscillation with two detectors in Kamioka and Korea 
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This paper updates and improves the study of electron neutrino appearance in the framework of two far 
detectors at different oscillation maxima, specifically, Tokai-To-Kamioka-to-Korea. We used a likelihood based 
on reconstructed quantities to distinguish charged current interactions from neutral current background. We 
studied the efficiency of the likelihood for a 20% photo-coverage in comparison of a 40% photo-coverage. We 
used a detailed neutrino event simulation to estimate the neutral current background. With these analysis tools 
we studied the sensitivity of the proposed experiment to CP violation and mass hierarchy as a function of the 
off-axis angle. 

PACS numbers: 14.60.Pq, 11.30.Er, 13.15.+g 



I. INTRODUCTION 

The 9i3 angle, the hierarchy of the largest mass splitting, 
and the CP violating phase are remaining undetermined quan- 
tities in the PMNS matrix formulation of neutrino masses and 
mixing ||Tl[2l- CP violation in the lepton sector and its re- 
lation to the matter anti-matter asymmetry of the universe is 
especially interesting. To determine the CP violating phase 
5 and the neutrino mass hierarchy, a powerful approach is to 
measure electron neutrino appearance at both the first and sec- 
ond oscillation maximum in a long baseline neutrino beam. 
Two different approaches have been considered in order to 
make this comparative measurement. One approach is to have 
two detectors at two different baselines of the same beam, 
one positioned for optimum response at the first oscillation 
maximum and the other positioned for optimum response at 
the second oscillation maximum. The optimum response is 
achieved using the off-axis technique |3] to produce a narrow 
energy band. This is the approach of the Tokai to Kamioka 
to Korea approach [4], henceforth referred to by the unofficial 
acronym T2KK. Another approach is to use an on-axis wide- 
band beam, and measure electron neutrino appearance from 
both the first and second maxima with a single detector fs*]. 
This is the approach employed by the BNL-FNAL working 
group [6] as a model for a long baseline neutrino oscillation 
experiment in the United States. 

In the first published T2KK article H, the off-axis angle of 
the Korean detector was assumed to be fixed at 2.5°. In this 
article, we study the sensitivity to CP violation and mass hi- 
erarchy if we choose a smaller off-axis angle for the location 
of the Korean detector, which blends the two approaches de- 
scribed above. As one can see in Fig.IT] the off-axis angle of 
1 .0° results in a fairly wide band beam, and we anticipate see- 
ing electron neutrino appearance at both the first and second 
maximum in the Korean detector The detector at the Kamioka 
location would remain at 2.5° off-axis, and be mainly sensi- 
tive to the first oscillation maximum. 

For this study, we assume an upgraded 1 .66 MW J-PARC 
beam created from 40 GeV protons, running 10^ seconds 
per year This is equivalent to 2.59 x 10^' protons-on-target 



(POT) per year. We assume 5 years of neutrino running and 
5 years of anti-neutrino running. The flux observed at four 
different off-axis angles, at 1050 km from the target, is shown 
in Fig. [T] We also assume two 0.27 Mton (fiducial volume) 
water Cherenkov detectors with 40% photo-coverage. One of 
them would be located at Kamioka, at a baseline of 295 km 
and at 2.5° off-axis angle from the beam. The second de- 
tector would be located in Korea at distances ranging from 
1000 to 1200 kilometers and off-axis angles ranging from 1° 
to 2.5°. The beam intensity assumed for this article is a factor 
of 2.4 lower than the beam assumed previously |4]. This more 
conservative beam power is being considered for benchmark 
studies 1 7]. In addition we assumed that a year of running is 
10^ seconds instead of 1.12 x 10^, so combining the change 
in beam intensity and running time, the number of POT per 
year is a factor of 2.7 lower than in previous studies. 

Our tool for these studies is the fully reconstructed at- 
mospheric neutrino Monte Carlo sample from the Super- 
Kamiokande experiment whereas the earlier paper used 
event rates calculated for the T2K experiment scaled for dis- 
tance. In order to simulate the T2KK beam, we re-weight 
events in the Super-K atmospheric neutrino Monte Carlo sam- 
ple by the ratio of the T2KK flux to the atmospheric flux. 



II. SIGNAL VS. BACKGROUND LIKELIHOOD ANALYSIS 



Our objective is to identify and reconstruct an excess of 
charged current Vg interactions in a nearly pure beam. We 
shall be especially interested in quasi-elastic interactions such 
as Vf.n e^p. In the experiment considered, the appearance 
probability is a few percent at most, and only a small number 
of events are anticipated above a non-negligible background. 
There are three categories of background: 

• Ve beam: The irreducible background from electron 
neutrinos in the beam flux regardless of neutrino oscil- 
lation. These come mainly from muon decay and Kej,. 

• Neutral current (NC): Background where the hadronic 
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FIG. 1 : (Color online) Neutrino flux as a function of energy for sev- 
eral off-axis angle, and a 0.75MW beam at 1050km from the target. 
For comparison, the probability, for two baselines consid- 

ered for T2KK (295km and 1050km). Neutrino mixing parameters 
are: normal hierarchy, Am^^i 31) ~ 8-0 x 10^^,2.5 x lO^-'eV^, and 

sin^ 26(12,23) = 0.86, 1.0. We take the earth density to be constant 
and equal to 2.8 g/cm?. 

recoil of neutral current interactions are misidentified as 
electron showers. 

• mis-ID: Background due to muons mis-identified as 
electron showers. 

The irreducible beam background is estimated from the de- 
tails of hadron production and muon decay in the beam Monte 
Carlo. We take as input the calculated flux from a beam simu- 
lation assuming a graphite cylinder target, 30 mm in diameter 
and 900 mm in length and a 130 m decay tunnel. These are 
the parameters of the T2K experiment. The neutral current 
background mainly consists of hadronic recoils with a single 
Tt*'. The Tt" decays into two photons and if one of the photons 
is missed because of a very small energy or an overlapping 
ring, then the event can be misidentified as a single electro- 
magnetic shower and therefore fake a Vg CCQE event. The 
dominant case is when one of the photons was missed be- 
cause the energy was too small. The v^, mis-ID background 
consists of charge current events where the Cherenkov ring 
from the outgoing muon is mis-identified as an electron by 
the reconstruction algorithm. This is the smallest source of 
background. 

Since we are interested in Vg appearance and especially 
undergoing quasi-elastic interactions, the events that we want 
to select are fully contained inside the fiducial volume, have a 
single Cherenkov ring identified as electron-like, and with no 
muon decay electron present (a muon decay electron would 



signal a missed 7t+ in a non-quasielastic final state). These 
are referred to as pre-cuts. Before building the likelihood, we 
applied these pre-cuts, in order to remove a significant part of 
the background. 

The pre-cut efficiencies are listed in Table |T] The NC effi- 
ciency is based on the total cross section for neutral current 
interactions which includes a large component of neutrino- 
nucleon elastic scattering. These are mostly unobserved in a 
water Cherenkov detector The NC background events that 
pass the pre-cuts are mostly single-Tt" production. 





Signal 


Background 


True V energy 


Ve (avg) 


QEv, 


non-QE 


NC 


mis-ID 


- 0.35 GeV 


95% 


94% 


53% 


0.4% 


0.5% 


0.35 - 0.85 GeV 


87% 


96% 


49% 


3% 


0.4% 


0.85 GeV - 1.5 GeV 


70% 


95% 


43% 


8% 


0.3% 


1.5 -2.0 GeV 


58% 


91% 


38% 


11% 


0.5% 


2.0 -3.0 GeV 


51% 


91% 


35% 


11% 


0.8% 


3.0 - 4.0 GeV 


45% 


90% 


34% 


12% 


0.9% 


4.0 - 5.0 GeV 


43% 


90% 


33% 


13% 


1.0% 


5.0 - 10.0 GeV 


37% 


86% 


29% 


10% 


1.4% 



TABLE I: Efficiency of pre-cuts as applied to neutrino interactions 
in the fiducial volume of the Super-Kamiokande detector simulation. 
The charged current Ve interactions are broken down separately for 
quasi-elastic and non-quasi-elastic samples. The NC sample includes 
elastic scattering in the denominator of the efficiency calculation. 



After applying pre-cuts, we make the final event selection 
using a likelihood based on several event characteristics and 
using the ROOT package TMVA J^. This is a similar ap- 
proach to one previously studied by others |10]. We recon- 
struct the neutrino energy assuming quasi-elastic interactions. 
This depends on particle masses, the reconstructed momen- 
tum and energy of the outgoing lepton, and the angle between 
the outgoing lepton direction and the known neutrino beam 
direction (9ve): 

mnEe~ml/2 

m„ —Ef.+ [Pe COS Svej 

The variables that are used in the likelihood can be divided 
into three categories: 

• Basic Super-Kamiokande event parameters: 

- The ring-finding parameter used to count rings 

- The e-like//j-like particle identification parameter 

• Light-pattern parameters used for vP finding: 

- The Tt*^ mass 

- The Tt*^ likelihood 

- The energy fraction of the 2"'' ring 

• Beam related variable: 
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- The angle between the outgoing lepton and the 
beam dkection 

We akeady cut on the ring parameter and the PID parameter in 
the set of pre-cuts (Table|l|. Here we used the continuous dis- 
tribution of these parameters as input to the Hkehhood. There 
are three variables related to a specialized fitter (POLf it for 
Pattern-Of-Light fitter) used to select single Tt'' events ifTTIl . 
The output of this fitter includes an overall likelihood as well 
as the best fit mass and energy fraction of the two gammas 
from T\P decay. We also use one variable that requires knowl- 
edge of the beam direction, and therefore is not a standard SK 
variable for atmospheric neutrino analysis. For that variable, 
we had to use the MC truth information about the neutrino 
direction in the simulated atmospheric neutrino Monte Carlo 
sample. Unlike the accelerator-based experiment, these events 
are simulated over a wide-range of incident angles. However, 
the Super-K detector has uniform response. The distributions 
of the combined likelihood for each energy bin is shown in 
FigE] The separation between signal and background is strik- 
ing at low energies but becomes worse at higher energies. 
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FIG. 2: (Color online) Combined likelihood distribution from 6 input 
variables, shown separately for 8 energy bins. Charged current Vg 
signal is shown in blue (filled), and the background is red (hatched). 
The events used have passed the defined pre-cuts. 

To choose where to cut on the likelihood variable, we com- 
pute the signal over square root of background, S/\/B for 
several positions of the cut. We tested cuts that range from 
keeping 10% of signal to keeping 100% of the signal (at the 
expense of increasing background). We also varied the off- 
axis angle and considered separate energy bins. We found 
that keeping a large fraction of signal, 80%, maximizes 5/\/B- 
The energy dependent efficiencies for an 80% likelihood cut 
is given in Table Ull 



A. Photo coverage 

Due to the accident that happened in November 2001, 
where about half of the Super-Kamiokande phototubes were 





Cut that keeps 80% of signal 


Energy (rec) 




NC 


mis-ID 


- 350 MeV 


80% 


15% 


15% 


350 - 850 MeV 


80% 


25% 


40% 


850 MeV - 1.5 GeV 


80% 


28% 


30% 


1.5 -2.0GeV 


80% 


30% 


32% 


2.0 - 3.0 GeV 


80% 


40% 


18% 


3.0 - 4.0 GeV 


80% 


50% 


28% 


4.0 - 5.0 GeV 


80% 


65% 


55% 


5.0 - 10.0 GeV 


80% 


45% 


18% 



TABLE II: Efficiency for the likelihood cut that keeps 80% of the 
signal. These efficiencies are calculated for events which have al- 
ready passed the pre-cuts, and are calculated based on reconstructed 
energy. 





\ 


e 


NC 


V 






Photo-cov. 


Photo-cov. 


Photo-cov. 


Energy (true) 


40% 


20% 


40% 


20% 


40% 


20% 


- 350 MeV 


95% 


93% 


0.4% 


0.3% 


0.5% 


0.5% 


350 - 850 MeV 


87% 


86% 


3% 


3% 


0.4% 


0.5% 


850 MeV - 1.5 GeV 


70% 


70% 


8% 


9% 


0.3% 


0.5% 


1.5 -2.0 GeV 


58% 


57% 


11% 


12% 


0.5% 


0.5% 


2.0 - 3.0 GeV 


51% 


49% 


11% 


12% 


0.8% 


0.6% 


3.0 - 4.0 GeV 


45% 


45% 


12% 


12% 


0.9% 


0.9% 


4.0 - 5.0 GeV 


43% 


45% 


13% 


13% 


1.0% 


0.9% 


5.0 - 10.0 GeV 


37% 


37% 


10% 


12% 


1.4% 


1.0% 



TABLE III: Pre-cut efficiency for two photo-coverage: 40% (SK-I) 
and 20% (SK-II). The NC sample includes elastic scattering in the 
denominator of the efficiency calculation. 



destroyed, Super-K has run with both 40% and 20% photo- 
coverage, and has atmospheric neutrino Monte Carlo samples 
for both conditions. It was therefore easy to repeat our studies 
for a detector with a 20% photo-coverage. Overall we found 
that both the pre-cuts and the likelihood are nearly as efficient 
in a detector with 20% coverage as they are in a detector with 
40% coverage. The comparison of the pre-cuts is presented in 
Table Uni and the comparison of the likelihood is presented in 
TableHYl 



III. HOW TO COMPUTE THE BACKGROUND 
SPECTRUM 

As mentioned in Section [III there are three categories of 
background: Vg beam background (Vg beam), neutral cur- 
rent background (NC), and charged current mis-identified 
background (v^ mis-ID). To simulate the background in the 
long baseline beam experiment, we used the SK atmospheric 
Monte Carlo as follows: 

• We ran over the atmospheric SK Monte Carlo, and kept 
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NC 




Photo-coverage 


Photo-coverage 


Energy (rec) 


40% 


20% 


40% 


20% 


- 350 MeV 


80% 


80% 


15% 


15% 


350 - 850 MeV 


80% 


80% 


25% 


24% 


850 MeV - 1.5GeV 


80% 


80% 


28% 


25% 


1.5-2.0GeV 


80% 


80% 


30% 


35% 


2.0 - 3.0 GeV 


80% 


80% 


40% 


40% 


3.0 - 4.0 GeV 


80% 


80% 


50% 


42% 


4.0 - 5.0 GeV 


80% 


80% 


65% 


50% 


5.0 - 10.0 GeV 


80% 


80% 


45% 


45% 



TABLE IV: Likelihood efficiency for two photo-coverages: 40% 
(SK-I) and 20% (SK-II). The likelihood cut keeps 80% of signal. 



IV. OFF-AXIS ANGLE ANALYSIS 

Using the cut on the hkelihood that keeps 80% of the signal, 
we present in Fig.|4]spectra at the Kamioka location and at the 
Korean location for 1° off-axis angle and 2.5° off-axis angle. 
We also present the sensitivity to mass hierarchy and CP vio- 
lation, for four different values of the off-axis angle position 
of the Korean detector. The analysis used to compute the 
sensitivity is similar to that previously used |3l and is defined 
as: 



k=\ V '=1 '^i 

where 
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FIG. 3: Smearing matrix for neutral current events: the result of 
energy reconstructed using Eqn.[T] versus true neutrino energy. 



events which passed all the pre-cuts. 

We applied the likelihood efficiency coiTesponding to 
the right background type (v,., mis-ID or NC) and 
using the reconstructed energy. This takes care of the 
likelihood efficiency, and also the energy resolution of 
the detector since we use reconstructed energy. 



• We re-weighted this background spectrum by the ratio 
of the beam flux to the atmospheric flux. 

• We normalized the final background spectrum in order 
to account for the running conditions of the experiment: 
volume of detector, beam power, etc. 

It is important to consider the neutral current background 
properly since its energy response is very uncorrected, as can 
be seen in Fig. [3] 



Here, Nexj) is the number of "experiments". For example if we 
have two detectors (Kamioka and Korea) and run with only 
neutrinos then A^^^ =2. If we have two detectors but run with 
neutrinos and anti-neutrinos then Nexp = 4. Compared to the 
publication of Ishitsuka et al, we added two energy bins and 
use events up to 3 GeV, which is relevant when the Korean de- 
tector is located at small off-axis angles. So for this analysis, 
we have Nexp = 4 since we ran for neutrinos and anti-neutrinos 
and have two detectors. We have 7 energy bins (NEhm)' 400- 
500 MeV, 500-600 MeV, 600-700 MeV, 700-800 MeV, 800- 
1200 MeV, 1200-2000 MeV, 2000-3000 MeV. The sum over 
j in Eq.|2]is the sum over the systematic errors. We consider 
fifteen systematic errors Gj in this study. They are presented 
in Table. [Viand they are split into three groups. The first seven 
errors are on the background, the next six on the signal, and 
the last two on the energy scale. The largest systematic un- 
certainty comes from the signal normalization above 1 .2 GeV, 
and this is due to the uncertainty on the number of rings for 
Multi-GeV electron-like events |12]. The systematic uncer- 
tainties were estimated using work by the Super-Kamiokande 
collaboration |8, 12, 13]. 

The results for the mass hierarchy and CP violation sensi- 
tivity are presented in Fig.|5]and Fig.|6] We find that the best 
sensitivity to both CP violation and mass hierarchy is achieved 
with the Korean detector located at 1° off-axis. The improve- 
ment in sensitivity to CP violation is rather minimal, however 
the sensitivity to mass hierarchy is improved by a factor of 
three compared to the original configuration with the Korean 
detector located at 2.5° off axis. This is due to the informa- 
tion gained by including the first oscillation maximum, with 
higher energy neutrinos, in the Korean far detector 

We note that several improvements have been made since 
the T2KK article published in 2005 [4]. Several minor prob- 
lems were fixed and the cut on the likelihood variable was 
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FIG. 4: (Color online) Reconstructed energy spectra at Kamioka (top), Korea 1.0° olf-axis (bottom left) and Korea 2.5° olf-axis (bottom 
right) for sin^(26i3) =0.04 and normal hierarchy. The remaining oscillation parameters are: Ara^jj = 8.0 X 10^^2.5 x lO^^eV^ and 

sin^ 26(12 23) = 0.86, 1.0. Each plot is normalized to 5 years of running with neutrino, a 1.66 MW beam with 40 GeV protons and in a 0.27 
Mton (FV) detector (i.e. 5 x 2.59 x 10^' POT). 



added. This allowed us to gain a significant number of signal 
events. For example in the 350-850 MeV bin, the combined 
efficiency (pre-cuts and likelihood) is 68%, where in the same 
bin of Ref. |4i] it was 40%. In addition, the likelihood cut al- 
lows us to increase S/Vb. Again for the 350-850 MeV bin, 
the S/Vb was increased by about 20%. If we had run with the 
same number of protons-on-target as the authors of the 2005 



paper, the sensitivity would be a factor of two better than what 
we are reporting for 2.5° off-axis angle. Conversely, with the 
conservative benchmark beam power of 1 .66 MW instead of 
4MW, our sensitivity with the Korean detector located at 2.5° 
off-axis is roughly equivalent to that of the 2005 paper 
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FIG. 5; Sensitivity to CP violation (left) and mass hierarchy (right) for different values of the off-axis angle. Other mixing parameters are the 
same as used in previous figures. Each plot considers 5 years of running with neutrinos and 5 years with anti-neutrinos, a 1.66 MW beam with 
40 GeV protons and in a 0.27 Mton (FV) detector, i.e. 10 x 2.59 x 10^' POT. The parameter region above the curve is the region where we 
can determine CP violation or mass hierarchy to either 2 or 3 O. 




FIG. 6: Sensitivity to CP violation (left) and mass hierarchy (right) for different values of the off-axis angle. Other mixing parameters and the 
beam exposure are the same as used in previous figures. These plots show for what fraction of possible values of the CP phase 5 we will be 
able to determine CP violation or mass hierarchy. The plain lines are for normal hierarchy while the dotted lines are for inverted hierarchy. 
If 013 is large enough the mass hierarchy can always be determined whatever the value of §; this is not the case for establishing CP violation 
when 5 approaches 0, it, or 2ir. 
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Index 


Systematic uncertainty 


Value 


1 


BG normalization below 1.2 GeV (for Kamioka) 


5% 


2 


BG normalization above 1.2 GeV (for Kamioka) 


5% 


3 


BG normalization below 1.2 GeV (for Korea) 


5% 


4 


BG normalization above 1.2 GeV (for Korea) 


5% 


5 


BG normalization between v^. and v^, below 1.2 GeV 


5% 


6 


BG normalization between v^. and v^, above 1.2 GeV 


5% 


7 


BG spectrum (common for Kamioka and Korea) 


5% 


8 


Signal normalization below 1.2 GeV o(v^)/ct(V(.) 


5% 


9 


Signal normalization above 1.2 GeV CT(v^)/o(Vf.) 


20% 


10 


[o(v^)/o(v,)]/[o(v7,)/o(v,)] below 1.2 GeV 


5% 


11 


[o(v^)/o(v,)]/[o(v7,)/o(v,)] above 1.2 GeV 


5% 


12 


Efficiency difference between Kamioka and 






Korea detector below 1.2 GeV 


1% 


13 


Efficiency difference between Kamioka and 






Korea detector above 1.2 GeV 


1% 


14 


Energy scale difference between Kamioka and 






Korea detector 


1% 


15 


Energy scale difference between near and 






(Kamioka/Korea) detector 


1% 



TABLE V: List of systematic uncertainties and their assumed values. 



V. CONCLUSIONS 

We have presented an updated and improved study of long 
baseline neutrino oscillation with a detector in Kamioka and a 



second in Korea. Using pre-cuts and a likelihood designed to 
reject neutral cuiTent background while keeping charged cur- 
rent quasi-elastic events, we were able to increase the amount 
of signal that we keep in the main signal bin (350-850 MeV) 
from about 40% to 68%, and we were able to remove more 
background than what was done before. We found that the 
effectiveness of the cuts and likelihood was relatively undi- 
minished when applied to a detector with 20% rather than 
40% photo-coverage. We found that the best location among 
the possibilities we explored for the Korean detector is 1.0°, 
which is more on-axis than previously considered, and allows 
a somewhat wider band neutrino energy spectrum. This im- 
proved the T2KK sensitivity by a factor of two compared to 
what was published previously, even after taking a more con- 
servative number of POT per year to be a factor of three lower. 
With an experiment configured at 1 .0°, and a benchmark beam 
power of 1 .66 MW, the neutrino mass hierarchy should be re- 
vealed if sin^29i3 is larger than 10^^ for a wide range of 5. 
CP violation would be detected at 3o significance for 70% of 
possible values of 8. 
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